Supernovae constraints on models of dark energy reexamined 
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We use the Type la Supernova gold sample data of Riess et al in order to constrain three models 
of dark energy. We study the Cardassian model, the Dvali- Turner gravity modified model and the 
generalized Chaplygin gas model of dark energy - dark matter unification. In our best fit analysis 
for these three dark energy proposals we consider flat model and the non-flat model priors. We 
also discuss the degeneracy of the models with the XCDM model through the computation of the 
so-called jerk parameter. 
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I. INTRODUCTION 



The surprising discovery of the present late-time accel- 
eration of the Universe [l| and the related fact that most 
of its energy is in the form of a mysterious dark energy is 
possibly one of the most puzzling issues of modern cos- 
mology. Several scenarios have been put forward as a 
possible explanation. A positive cosmological constant, 
although the simplest candidate, is not particularly at- 
tractive given the extreme fine tuning that is required 
to account for the observed accelerated expansion. This 
fact has led to models where the dark energy component 
varies with time, such as quintessence models ^. In these 
models, the required negative pressure is achieved trough 
the dynamics of a single (light) scalar field fs^ or, in some 
cases, two coupled scalar fields 4]. Despite some pleasing 
features, these models are not entirely satisfactory, since 
in order to achieve fix ~ (where fix and flm are 
the dark energy and matter energy densities at present, 
respectively) some fine tuning is required. Many other 
possibilities have been considered for the origin of this 
dark energy component such as a scalar field with a non- 
standard kinetic term and k-essence models ^5] ; it is also 
possible to construct models which have wx = p/ P < —1, 
the so-called phantom energy models Q. 

Recently, it has been proposed that the evidence for a 
dark energy component might be explained by a change 
in the equation of state of the background fluid, with 
an exotic equation of state, the generalized Chaplygin 
gas (GCG) model _7, ,8, 9]. The striking feature of this 
model is that it allows for an unification of dark energy 
and dark matter 

Another possible explanation for the accelerated ex- 
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pansion of the Universe could be the infrared modifi- 
cation of gravity one should expect from, for instance, 
extra dimensional physics, which would lead to a modifi- 
cation of the effective Friedmann equation at late times 
[Tol IT^ . An interesting variation of this proposal 
has been suggested by Dvali and Turner lOj (hereafter 
referred to as DT model). Another possibility, also origi- 
nally motivated by extra dimensions physics, is the mod- 
ification of the Friedmann equation by the introduction 
of an additional nonlinear term proportional to p", the 
so-called Cardassian model • 

Currently type la supernovae (SNe la) observations 
provide the most direct way to probe the dark energy 
component at low to medium redshifts. This is due to 
the fact that supernova data allows for a direct measurent 
of the luminosity distance, which is directly related to the 
expansion law of the Universe which, in turn, is the phys- 
ical quantity that is directly related with dark energy or 
that is modified by extra dimensional physics. This ap- 
proach has been explored by various groups in order to 
obtain insight into the nature of dark energy. Indeed, 
recently supernova data with 194 data points has been 
analysed [lj| and it was shown that it yields relevant con- 
straints on some cosmological parameters. In particular, 
it is possible to conclude that, when one considers the 
full supernova data set, the decelerating model is ruled 
out with a significant confidence level. It is also shown 
that one can measure the current value of the dark en- 
ergy equation of state with higher accuracy and the data 
prefers the phantom kind of equation of state, wx < ~1- 
Furthermore, the most significant result of that analysis 
is that, without a flat prior, that supernova data does 
not favor a flat ACDM model at least up to 68% confi- 
dence level, which is consistent with other cosmological 
observations. In what concerns the equation of state of 
the dark energy component, it has been shown in Ref. 
[13 , using the same set of supernovae data, that the best 
fit equation of state of dark energy evolves rapidly from 
Wx — in the past to wx ^ ~1 in the present, which 
suggests that a time varying dark energy is better fitted 
with the data than the ACDM model. This result is also 
robust to changes of fim and remains valid for the inter- 
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val 0.1 < flm ^ 0.5. Supernova data has also been used 
in the context of different cosmological models for dark 
energy [H [13, El E 113, 113 ■ In this paper, we analyze 
the Cardassian, the DT and the GCG models in light 
of the Riess et al. SNe la compilation of data ,2^] . We 
consider both flat and non-flat priors. 

Notice however, that our analysis is restricted to the 
very late history of the Universe and does not address 
dark energy effects on the cosmic microwave background 
fluctuations or on structure formation. 

This paper is organized as follows. In section 11 we 
describe our best fit analysis of SNe la data that will be 
employed as methodology to constrain the dark energy 
models we have studied, namely the Cardassian, DT and 
GCG models. In section III we consider the best fit anal- 
ysis to constrain the Cardassiam dark energy model. In 
section IV, we discusss the DT model and in section V 
the dark energy - dark matter unification GCG model. 
Section VI is devoted to the discussion of the degeneracy 
of the discussed models with the XCDM model through 
the introduction of the so-called jerk parameter. In sec- 
tion VII we present our conclusions. 



Parameters 


Range 


r^rn or As 


10, i[ 


nk 


]-!,![ 


M' 


[41,45] 


Cardassian: n 


[-30,2/3] 


DT: /? 


[-60, 1] 


GCG: a 


[0, 10] 



TABLE I: Parameter range for the best-fit analysis 



was not certain or the photometry was incomplete, in- 
creasing the reliability of the sample. The second set 
contains the 143 points from the first one plus 14 points 
discovered recently using Hubble Space Telescope (HST) 
[2^ and is named as gold sample, following Riess et al. 
We will name the first sample as gold w/o HST. The 
main difference between the two samples is that the full 
gold sample also covers higher redshifts (1 < z < 1.6). 

The data points in these samples are given in terms of 
the distance modulus 



II. OBSERVATIONAL CONSTRAINTS FROM 
SUPERNOVAE DATA 

The observations of supernovae measure essentially the 
apparent magnitude m, which is related to the luminosity 
distance (ii by 



fJ-ohs{z) = m{z) - Mohs{z) , 



(6) 



where 



m{z) = M + 5 logio Dl (z) 



Dl{z) = ^dL{z) 
c 



(1) 



(2) 



is the dimensionless luminosity distance and 

dL{z) = (1 + z)dM{z) , (3) 
with d},[{z) being the comoving distance given by 

1 



Also, 



duiz) 



M + ^ log 



H{z') 



-dz' 



10 



1 Mpc 



25 



(4) 



(5) 



where M is the absolute magnitude which is believed to 
be constant for all supernovae of type la. 

For our analysis, we consider the two sets of supernovae 
data recently compiled by Riess et al. |22] ■ The first set 
contains 143 points from previously published data that 
were taken from the 230 Tonry et al. 23] data alongwith 
the 23 points from Harris et al. |2^ . They have discarded 
various points where the classification of the supernovae 



and the errors fJ^^-i^Sz) already quoted take into account 
the effects of peculiar motions. 
The is calculated from 



E 



Mobs(2i) - TW' - 5 logio DLth{zi; Ca) 



, (7) 



where M' 



Ai — Mobs is a free parameter and 



DLth{z', Ca) is the theoretical prediction for the dimen- 
sionless luminosity distance of a supernova at a partic- 
ular distance, for a given model with parameters Cq. It 
can be computed for each model from the Friedmann ex- 
pansion law (cf. Eqs. ©, lfn|) or lfT7|) below) combined 
with Eqs. 

In the following, we will consider the three models re- 
ferred to in the Introduction and perform a best fit analy- 
sis with the minimization of the x^, Eq. |(7J), with respect 
to M' , Qrm and the respective model parameter (s), 
using a MINUIT ^ based code. 

The allowed variation range of the parameters is pre- 
sented in Table HJ Al' is a model independent parameter 
and hence its best fit value should not depend on the spe- 
cific model. We found that the best fit value for A4' for 
all the models considered here is 43.3, which is consistent 
with the one obtained in [ij. Hence, we have also used 
this value for Ai' throughout our analysis. We have also 
checked that the result does not change if we marginalise 
over Ai'. 

In what follows we shall present a description of the 
models that are considered in this analysis and perform 
the best fit study considering flat and non-flat priors. 
Our results are summarized in Table HTl 
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Data Sample 


^rn or As 


Model parameter 






Cardassian model 


Flat 


Gold w/o HST 


0.53 


-2.0 




154.6 


Prior 


Gold 


0.49 


-1.4 




173.7 


Non-Flat 


Gold w/o HST 


0.97 


-0.93 


-0.75 


154.4 


Prior 


Gold 


0.21 


-3.1 


0.47 


173.2 


DT model 


Flat 


Gold w/o HST 


0.55 


-60.0 




155.4 


Prior 


Gold 


0.51 


-19.2 




174.7 


Non-Flat 


Gold w/o HST 


1.0 


-15.6 


-0.71 


155.0 


Prior 


Gold 


0.24 


-60.0 


0.43 


174.0 


GCG model 


Flat 


Gold w/o HST 


0.98 


6.2 




155.0 


Prior 


Gold 


0.93 


2.8 




174.2 


Non-Flat 


Gold w/o HST 


0.73 


1.3 


-1.0 


154.2 


Prior 


Gold 


0.97 


4.0 


0.02 


174.5 



TABLE H: Best fit parameters for the Cardassian, DT and GCG models, for the two data samples of Riess et al, considering 
flat and non-flat priors. The best fit value used for A4' is 43.3. 




FIG. 1: Confidence contours in the Qm — n parameter space for the fiat Cardassian model. The solid and dashed lines represent 
the 68% and 95% confidence regions, respectively; in the left panel are shown the results for the gold sample without the HST 
SNe la, whereas in the right panel, the full gold sample is taken into account. 



III. CARDASSIAN MODEL 

We first consider the so Cardassian model which 
justifies the late time accelerating universe by invoking an 
additional term in the Friedmann equation proportional 
to p". In this model, the universe is composed only of 
radiation and matter (including baryon and cold dark 
matter) and the increasing expansion rate is given by 

where Mpi = 1.22 x lO^^GeV is the 4-dimensional Planck 
mass, b and n are constants, and we have added a curva- 



ture term to the original Cardassian model. At present, 
the universe is matter dominated, i.e. pm ^ Prad, hence 
P ~ Pm ■ The new term dominates only recently, at about 
z ^ 1, hence in order to get the recent acceleration in the 
expansion rate, n < 2/3 is required. 

The theoretical motivation for the Cardassian term 
is fairly speculative. It can be argued that its origin 
may arise as a consequence of embedding our (3+1)- 
dimensional brane universe in extra-dimensions |26j| . or 
from some unknown interactions between matter parti- 
cles 113 . 

In a matter dominated universe, equation (|HJ| can be 



4 




FIG. 2: Confidence contours in tlie and Q,n — n, Qk ^ and Q.^ — n parameter space for tlie non-flat Cardassian model. As 
in Figure Q tiie solid and dashed lines represent the 68% and 95% confidence regions, respectively; in the left panel is used the 
golden sample without HST SNe la, whereas in the right one the full gold sample is taken into account. 



rewritten as where Hq is the present day value of the Hubble constant 

and = — 3773- is the present curvature parameter. 



H 



-- rir)i(l + z)^ + f^fe(l + z)"^ Notice that the only parameter of the model is n and 

-(l-a„-f]fc)(l + z)3" , (9) 
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FIG. 3: Confidence contours in the fim — P parameter space for the fiat DT model. The solid and dashed lines represent the 
68% and 95% confidence regions, respectively; in the left panel are shown the results for the gold sample without the HST SNe 
la, whereas in the right the full gold sample is taken into account. 



that the case n — Q corresponds to the ACDM model. 

The best fit results for this model, for the two data 
samples we are considering, and taking into account flat 
and non-flat priors, are summarized in Table HTl 

For the flat case, we have only two parameters and the 
best fitting results we get are {ri„i,n} — {0.53,-2.0}, 
without HST, and {n^,n} = {0.49,-1.4}, with HST 
SNe la. In Fig. ^we show the 68% and 95% confidence 
contour plots. These results are consistent with those 
obtained in previous works by Gong & Duan ^sf , Zhu et 
al. and Sen & Sen for other data sets. Hence, 
we see that the case n = 0, which corresponds to the 
ACDM model, is excluded at a 95% confidence level by 
both samples, even though HST sample favors slightly 
larger values of n. 

If we relax the flat prior and consider the curvature we 
find that the best fit analysis reveals that the gold sample 
without HST data favors a negative curvature around 
fife — —0.75, with a significant matter component, — 
0.97, even though the 68% confidence contour region (see 
Fig. 121) is consistent with curvatures in the range [—1,0.1] 
and matter density in the range [0.38,1]; for the model 
parameter we find that it lies in the range [—2.4, —0.2] 
with central value n = —0.93. From the 95% confidence 
contour, meaningful bounds can be still obtained for the 
model parameter: n must lie in the range [—3.6,0.1], 
hence not excluding the ACDM model. Values for Q.^ 
and, mainly, cannot be significantly constrained. 

Clearly the HST data brings the amount of matter to 
lower values, but pushes curvature to positive values and 
the model parameters for values that are smaller than 
the ones obtained with gold without HST data. We find 
that the best fitting value for the curvature is signifi- 
cantly positive Ofc — 0.47 and = 0.21; for the model 
parameter we obtain as best fit value —3.1. Moreover, 
the n = case is excluded with 95% confidence level. 



The contour plots are presented in Fig|2 



IV. DVALI-TURNER MODEL 

The second model we will consider is the one proposed 
by Dvali & Turner [lol ]. where the Friedmann equation 
is modified by the addition of the term /r^~^ , which 
can arise in theories with extra dimensions |ll|: Tc is a 
crossover scale which sets the scale beyond which the 
laws of the 4-dimensional gravity breakdown and be- 
come 5-dimensional. In this case, in contrast with the- 
ories with infinite volume extra dimensions, the laws of 
gravity are modifled in the far infrared and the cosmo- 
logical evolution gets modifled at late times; the short 
distance gravitational dynamics is very close to that of 
the 4-dimensional Einstein gravity, hence the early times 
cosmological evolution is very close to the Friedmann- 
Robertson- Walker picture. 

As a motivation for this type of modification of gravity, 
consider the model with a single extra dimension pdlll^. 
with the effective low-energy action given by 

S = d^xdyy^TZ 
rc J 

+ J d^x^ (M^i R + Csm) , (10) 

where y is the extra spatial coordinate, g^^^ is the trace 
of the 5-dimensional metric g'^g (^, -S — 0, 1, 2, 4), 
g the trace of the 4-dimensional metric induced in the 

brane, gf^u{x) = gl^J{x,y = 0). The first term in 
the action is the bulk 5-dimensional Einstein action, 
where TZ is the 5-dimensional Ricci scalar, and the sec- 
ond one is the 4-dimensional Einstein term, localized on 
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FIG. 4: Confidence contours in the Qm — P and Q.k — Q.m parameter space for the non-flat DT model using the gold sample 
without HST SNe la (left panel) and in the Q.m — P and Qk — (3 plane using the full gold sample (right panel). As in Figure |3] 
the solid and dashed lines represent the 68% and 95% confidence regions, respectively. 



the brane (at y — 0), where R is the 4-dimensional Ricci 
scalar and Csm is the Lagrangian of the fields in the 
Standard Model. For the maximally symmetric FRW 
ansatz, dsl — f{y,H)ds\ — dy^, where ds^ is the 4- 
dimensional maximally-symmetric metric, one gets the 
modified Friedmann equation on the brane 

± — ^ , (11) 

where p is the total energy density in the brane. 

Inspired in this construction, Dvali & Turner consid- 
ered a more generic (and radical) modification of the 
Friedmann equation [lOj 

rr2 8n k 

where we have also introduced the curvature term in the 
brane. The crossover scale Tc is fixed in order to eliminate 



the need for dark energy, 

= H^' (1 - - , (13) 

where we assume again that the universe is matter domi- 
nated. The Friedmann expansion law can then be written 
as 

="™(i + ^)' + (i-"'»-^fc) (;^) 

+Qkil + zf. (14) 

Notice that /3 is the only parameter of the model: for 
/3 = the new term behaves like a cosmological con- 
stant, and for /3 = 2 it corresponds to a "renormaliza- 
tion" of the Friedmann equation. Note also that the 
case /3 — I corresponds to the model in Eq. ^2 here- 
after called Dvali-Gabadadze-Porrati (DGP) model [lll |. 
The successful predictions of the Big-Bang nucleosynthe- 
sis impose a limit on /3, namely, f3 < 1.95; a more strin- 
gent bound follows from requiring that the new term does 
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not interfere with the formation of large-scale structure: 
/3 < 1. Moreover, it can be shown that, in the recent past 
(10'* > z 3> 1), this correction behaves like dark energy 
with equation of state We// = — 1 + /3/2, and w = — 1 in 
the distant future and, moreover, it can mimic w < —\ 
without violating the weak-energy condition 

For the flat DT model, the best-fitting parameters are 
= {0.55,-60.0} without HST and = 
{0.51, —19.2} with HST supernovae, as is summarized in 
Table HU From Fig. El where the 68% and 95% confi- 
dence contour plots are shown, it is clear that (5 is very 
weakly constrained by both supernovae data sets, and 
can become arbitrarily large and negative. Moreover, the 
cosmological constant, corresponding to /3 = 0, seems to 
be disfavoured; also /? = 1, the DGP model, is strongly 
disfavoured. These results are consistent with those of 
Elgaroy et al. obtained with another data set. 

If we fix = 1 (DGP model), and allow only ri,„ to 
vary, the best fit value are flm — 0.17 and — 0.16 
for the gold sample and gold sample without HST SNe 
la, respectively, which are consistent with the results of 
Gong & Duan 18]. 

Now, if we relax the fiat prior, we find that {flk, (3} ~ 
{—0.71, —15.6} are the best fit values for the gold sample 
without HST SNe la (Table HJ); the best fit value for fJ™ 
is in the end of the variation range we considered for this 
parameter, Qm = 1- As we can see from the contour 
plots shown in Fig01 cannot be constrained at 95% 
confidence level, however the 68% confidence contours 
show that data prefer f2m > 0.25. In what concerns /3, 
it can become arbitrarily large and negative, if we allow 
rim to be large. Moreover, both ACDM and DGP models 
are disfavored with 68% confidence. 

Once more the HST data brings the amount of 
matter to lower values and pushes curvature to pos- 
itive values. The best fit results are {D,m,^k, P} — 
{0.24,0.43,-60,0}. Also in this case both ACDM and 
DGP models are disfavored. 



V. GENERALIZED CHAPLYGIN GAS MODEL 

Finally, we consider the generalized Chaplygin gas 
model, which is characterized by the equation of state 



Pch 



A 

Pch 



(15) 



where A and a are positive constants. For a ~ I, the 
equation of state is reduced to the Chaplygin gas scenario 



eq 



Integrating the energy conservation equation with the 
equation of state lfT5|l . one gets [Q] 



Pch — PchO 



A, 



a3(l+a) 



l/(l+a) 



(16) 



where pcho is the present energy density of GCG and 
As = A/pl^f^^K One of the most striking features of this 



expression is that the energy density of this GCG, pch, in- 
terpolates between a dust dominated phase, pch oc a^'^, 
in the past and a de-Sitter phase, pch — ~Pch, at late 
times. This property makes the GCG model an interest- 
ing candidate for the unification of dark matter and dark 
energy. Moreover, one can see from the above equation 
that As must lie in the range Q < As < 1: for As — 
GCG behaves always as matter whereas for = 1, it 
behaves always as a cosmological constant. We should 
point out, however, that if we want to unify dark matter 
and dark energy, one has to exclude these two possibil- 
ities resulting the range for As as Q < As < 1. Notice 
also that a = corresponds to the ACDM model. 

The Friedmann equation for a non-flat unified GCG 
model in general is given by 



— ={l-nk) A, + (l-A,)(l + z)3(i+") 

-Ho/ L J 

+nkii + zf . (17) 



This model has been thoroughly scrutinized from 
the observational point of view; indeed, its compatibil- 
ity with the Cosmic Microwave Back grou nd Radiation 
(CMBR)peak location and amplitudes [H 113 , with SNe 
la data 0,|1O| and gravitational lensing statistics [sills^ 
has been analysed b y va rious groups. The issue of struc- 
ture formation |^ l^lsp and its dificulties l3j| have been 
recently addressed More recent analysis, based on 

the 194 SNe la data points from Ref. has yielded in- 
teresting results in what concerns the allowed parameter 
space of the model 

1) Data favors a > 1, although there is a strong degen- 
eracy on a. At 68% confidence level the minimal allowed 
values for a and As are 0.78 and 0.778, which rules out 
the ACDM model a = case. However, at 95% confi- 
dence level it is found that there is no constraint on a. 

2) Dropping the assumption of fiat prior, it is found that 
GCG is consistent with data for values of a sufficiently 
different from zero. Allowing for some small curvature, 
positive or negative, one finds that the GCG model is a 
more suitable description than the ACDM model. 

These results are similar to the ones obtained in Refs. 
0, lisll. where it is concluded that the supernova data of 
Ref. '23] favors "phantom" -like matter with an equation 
of state of the form p = ujp with w < — 1. 

In our present analysis, we have new results both with 
and without fiat prior. For the flat case, the 68% and 95% 
confidence level contours are quite similar both with and 
without HST data with that obtained by ,17] as shown 
in Fig. But still there is one interesting feature to 
note. It appears now that without the HST data, ACDM 
model [a — 0) is ruled out even at 95% confidence level 
which was not the case in the previous analysis. But 
incorporating the HST data, again makes ACDM model 
consistent at 95% C.L although it is still ruled out at 68% 
C.L. 
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When we relax the condition of flat prior, one can have some interesting features (see Fig. 6). Without the HST 
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data, supernova data prefers a nonflat model, but flat 
case is consistent both at 68% and 95% C.L if one choose 
a value for a sufficiently different from zero. This means 
without flat prior, and without HST data, data prefers 
a curved model for GCG but flat case is still consistent 
which is similar to the conclusion in p^ . 

But now if one includes the HST data, the best flt 
model itself becomes a very close to flat case {ilk — 0.02, 
Table nj which is completely a new feature. This shows 
that with the new HST data, the current gold sample 
of supernova data prefers a flat GCG model which is 
consistent with CMBR observation. It also shows that 
GCG is a better choice among the three possible choices 
of alternative models that we have considered here as far 
the gold sample of supernova data including those from 
HST measurements are concerned. 



VI. DEGENERACY WITH XCDM MODEL 

To illustrate the degeneracy between our models and 
the XCDM model, with constant equation of state P/ p = 
wx and dark energy density fix = l~i^m, let us consider 
the Taylor expansion of the luminosity distance as 



dL = 



c 



z + 



(1 - go) 



(18) 



where go is the deceleration parameter, related with the 
second derivative of the scale-factor, and jo is the so- 
callled statefinder or jerk parameter [s^, related with 
the third derivative of the scale-factor. The subscript 
"0" means that quantities are evaluated at present. The 
jerk parameter is related with the deceleration parameter 
qo as 



jo = go + 2(7o + 



dq 

dz 



(19) 



Table 



Notice that the jerk paramter is related to the geometry 
of the Universe (see Ref. [s^, where it is shown that 
the measurement of the cubic correction to the Hubble 
law via high-redshift supernovae is the flrst cosmological 
measurement, after the CMBR, that probes directly the 
effects of spatial curvature) . We have calculated go and 

for our models and obtained the results shown in 

HII We have considered only the flat case and, 
for simplicity, we study the DGP model instead of the 
more generic DT model; 51^ is the only parameter of 
this model. 

For low redshifts, it is sufficient to consider the first two 
terms in the series expansion of the luminosity distance, 
Eq. l(TF|) . From the expression of go for the GCG, we can 
see that in this case the SNe la can only constrain As, 
as go is independent of a. However for the Cardassian 
model the go dependence is both on fl^a and n, allowing 
low redshifts SNe la to put constraints in the model. 



Moreover, in order to have degeneracy between the 
models we are analyzing and XCDM model, the go 
parameter of these models must be equal, q^'-^^^^ — 



Model i .jyi^ich results that: 



Wx ^ (n- 1) 



1 - SI 



Qx 

1 



wx = - 



Wx 



(1 + n„^)nx 
nx 



for Cardassian model ; (20) 
for DGP model ; (21) 
for GCG model . (22) 



If one goes to higher redshifts, one also has to consider 
the higher order terms in the expansion of ^^(z); hence, 
in this case also the jerk parameters have to be equal, 



which means that ^ 

dz 

for Cardassian model 



have to be the same. We get that 



wx 
fix 



n 
1 - 



- 1 

fir. 



(23) 
(24) 



meaning that the dynamical evolution of this model is 
equivalent to a dark energy model with the same matter 
density and a constant equation of state given by Wef / — 
n — 1; the equivalent dark energy potential can be written 
&sV{<j)) = A[smhk{<j)/a + C)]'" with a ^ -2-2/(n-l) 
[Tsj (see also Ref. ^). We have seen that negative 
values of n are preferred which is consistent with the fact 
that phantom e quat ion of state, wx < —1, is favored by 
the data [11 El 113 • 

For DGP model one finds 



fl^ 



2 jf-f 1 



Wx 



fix = 



[l + flraf 
fll - 1 



1 



(25) 
(26) 



We see that for flm < 0.6 the the equation of state is 
phantom-like and, as f2„j — > 0, wx — > —1- We saw that 
the SNe la prefer fl„i ^ 0.2, consistent with the fact 
that phantom equation of state is favored. Moreover, 
the amount of matter needed in a XCDM model to be 
equivalent to a given DGP model is larger, fl^^^^^^ = 
1 — fix > flm- 

For the GCG model we find 



Wx = ~a{l - As) ~l , 
fix = 



l + a(l-^) ' 



(27) 
(28) 



We see that for any GCG model, the corresponding dark 
energy rnodel equation of state has to be always phantom 
type [13, as a > and < As < 1. Nevertheless, for 
low redshift data, GCG is degenerate with all kinds of 
constant equation of state dark energy model, including 
ACDM model, as can be inferred from Eq. if^ . 
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Model 




90 


dz 


n 


XCDM 


5 + 






Cardassian 


i + |(n 


- i)(i 


|(n- l)2n„,(l-n„) 


DGP 






9f2n.(l- 




2 i+n™ 


(l+O 




GCG 


la- 




|A(l-^s)(l + a) 



TABLE III: Deceleration and jerk parameters for XCDM, Cardassian, DGP and GCG models. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have performed likehhood analysis 
of the latest type la supernova data for three distinct 
dark matter models. We have considered the Cardassian 
model, the modified gravity Dvali- Turner model and the 
generalized Chaplygin gas model of unification of dark 
energy and dark matter. We find that SNe la most re- 
cent data allows, in all cases, for non-trivial constraints 
on model parameters as summarized in Table IL We find 
that for all models relaxing the flatness condition implies 
that data favors then a considerable negative curvature 
for the gold without HST dataset. For the gold dataset 
the resulting best fit value for the curvature is positive 
(the GCG model is nearly flat in this case). For all mod- 
els we have found, in what concerns the deceleration and 



jerk parameters, the conditions under which they are de- 
generate to the XCDM model. Thus, SNe la data clearly 
favors phantom-like equivalent equations of state. 

Finally, in what concerns the gold sample of supernova 
data including those from HST measurements, our anal- 
ysis reveals that that the GCG flat model is the better 
choice among the three possible alternative models that 
we have considered in this paper. 
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